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Many fundamental one-dimensional lattice models such as the Heisenberg or the Hubbard 
model are integrable. For these microscopic models, parameters in the Luttinger liquid 
theory can often be fixed and paramctcr-frcc results at low energies for many physical 
quantities such as dynamical correlation functions obtained where exact results are still 
out of reach. Quantum integrable models thus provide an important testing ground for 
I low-energy Luttinger liquid physics. They are, furthermore, also very interesting in their 

own right and show, for example, peculiar transport and thcrmalization properties. The 
Ch consequences of the conservation laws leading to integrability for the structure of the 

low-energy effective theory have, however, not fully been explored yet. I will discuss the 
connection between integrability and Luttinger liquid theory here, using the anisotropic 
Heisenberg model as an example. In particular, I will review the methods which allow to 
I fix free parameters in the Luttinger model with the help of the Bethe ansatz solution. As 

^ applications, parameter-free results for the susceptibility in the presence of non-magnetic 

^v^j impurities, for spin transport, and for the spin-lattice relaxation rate are discussed. 
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1. Introduction 

The Tomonaga-Luttinger liquid 1-4 is believed to describe the low-energy properties 
t-H of gapless one-dimensional interacting electron systems irrespective of the precise 

nature of the microscopic Hamiltonian. This universality can be understood in a 
• renormalization group sense as irrelevance of band curvature and additional inter- 

action terms which might arise when deriving this low-energy effective theory from 
a microscopic model. Similar to the important role Onsager's exact solution 5 of 
the two-dimensional Ising model has played in establishing and confirming general 
renormalization group theory, exactly solvable one-dimensional quantum models 
have been crucial for the development of Luttinger liquid theory. 

Integrable models are, furthermore, also interesting in their own right and a 
number of almost ideal realizations are known today. One example are cuprate 
spin chains such as Sr 2 Cu0 3 whose magnetic properties are well described by the 
integrable one-dimensional Heisenberg model. 6-14 Furthermore, cold atomic gases 
represent quantum systems which are to a high degree isolated from the surround- 
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ings and whose Hamiltonians are easily tunable. This makes it possible to use them 
as quantum simulators to study almost perfect realizations of integrable systems 
such as the Lieb-Liniger model 15,16 or the fermionic Hubbard model. 17 

In Sec. 2 I will discuss quantum integrability with a particular emphasis on Bethe 
ansatz integrable models and outline possible effects on transport and the thermal- 
ization of closed quantum systems. In the rest of the paper, I will then concentrate 
on the anisotropic Hcisenbcrg (or XX Z) model as one specific example for a Bcthc 
ansatz integrable model. In Sec. 3 I describe how the Luttinger model, including 
leading irrelevant operators, can be obtained from this microscopic model by us- 
ing bosonization techniques. In Sec. 4 I then briefly outline important aspects of 
the Bethe ansatz solution. In Sec. 5 it is shown that a comparison of the results 
of Sec. 3 and Sec. 4 allows to fix parameters in the Luttinger liquid theory for 
the XX Z model. Applications of the parameter-free low-energy effective theory to 
calculate various properties of spin chains are considered in Sec. 6. This includes 
the calculation of susceptibilities in the presence of non-magnetic impurities, and 
results for spin transport and NMR relaxation rates. The final section is devoted to 
a brief summary and some conclusions. 



2. Quantum integrability 

A classical system with Hamilton function H and phase space dimension 2N is 
integrable if it has N constants of motion Q n with 

{H,Q„} = and {Q i ,Q j } = if i^j. (1) 

Here {., .} denotes the Poisson bracket. Quantum integrability, on the other hand, is 
much harder to define precisely, see, for example, Ref. 18. In this regard it is impor- 
tant to note that every quantum system in the thermodynamic limit, irrespective 
of integrability, has infinitely many conservation laws 

[H,\E n )(E n \}=0 (2) 

where H is the Hamiltonian with eigenstates \E n ) and [., .] denotes the commutator. 
Apart from these non-local conservation laws a quantum system can have local 
conservation laws given by 

Qn = ^2q n , 3 or Q n = dxq n (x) (3) 
j J 

where q n j is a density operator acting on n neighboring sites in the case of a lattice 
model while for a continuum model q n (x) is a fully local density operator. A generic 
example for a local conservation law is the Hamiltonian itself for models with short 
range interactions. In Bethe ansatz integrable models a whole set of such local 
conservation laws does exist which can be obtained from the transfer matrix of the 
corresponding two-dimensional classical model by taking successive derivatives of 
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the transfer matrix r with respect to the spectral parameter A 

(4) 

Here £ is the spectral parameter at which the transfer matrix is evaluated. These 
conserved quantities are directly related to the existence of so-called i?-matrices 
which fulfill the Yang-Baxter equations and from which the transfer matrices r(A) 
can be constructed. 19 

For a low-energy effective theory describing such an integrable model, we have 
to demand — at least in principle — that the low-energy Hamiltonian H e s also fulfills 

[fleff, Qn] = 0. (5) 

This corresponds to a fine-tuning of parameters in the Luttinger model. In particu- 
lar, it might mean that certain terms which are not forbidden by general symmetry 
considerations have to vanish. Such a program has not fully been explored yet; in 
Sec. 5.2 we will see, as an example, that the conserved quantity Q3 for the XXZ 
model does indeed prevent certain terms from occuring in the low-energy theory. 



Sn0C dF lnr(A) 



2.1. Consequences for transport 

Local conservation laws can have a dramatic effect on the transport properties. 20 
This can be easily understood as follows. We can always define a local current 
density ji by making use of the continuity equation 

d 

-Q t Pi + 3l - Ji-i = (6) 

where pi is the density at site I. The current itself is then given by J = J2i3i- This 
current could be, for example, an electric, spin or thermal current. Conserved quan- 
tities Q n can now prevent a current from decaying completely leading to ballistic 
transport and a finite Drude weight 

Here Q n can denote a local or non-local conserved quantity, L is the length of 
the system, and T the temperature. The second relation in Eq. (7) is the Mazur 
inequality 21 which becomes an equality if all conservation laws, local and non-local, 
arc included. 22,23 In order to obtain a possible non-zero Drude weight at finite 
temperatures within a Luttinger model description, the relevant conservation laws 
have to be taken into account explicitly. One way to achieve this is discussed in 
Sec. 6.2. Importantly, one expects that only local or pseudo-local conservation laws 
Q n with (JQ n ) ^ can give rise to a finite bound in Eq. (7) so that D(T > 0) 7^ 
is characteristic for an integrable model. 
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2.2. Consequences for thermalization in closed systems 

Additional local conservation laws can also have a profound impact on a possible 
thermalization of a closed quantum system. Imagine that we prepare an initial 
state |^(0)) and follow the unitary time evolution of this state under an integrable 
Hamiltonian. One says that a closed quantum system in the thermodynamic limit 
has thermalized if for any local observable O the limit 

Ooo= lim<*(t)|0|*(t)) (8) 

t— too 

is well-defined and time independent and can also be expressed as an ensemble 
average 

Ooc = Tr(Op) (9) 

with an appropriately chosen density matrix p. Note that even for a generic closed 
quantum system temperature T is not defined by an external bath but rather by 
the energy of the initial state 

(tf(0)|#|tf(0)) =Tr(Hc- H / T )/Z (10) 

v v ' 

Tt( Pc h) 

with T acting as a Lagrange multiplier and Z being the partition function. For an 
integrable model we have to demand that the relation (10) also holds if we replace 
H — > Q n and the canonical density matrix p c by the density matrix 24 

p=Iexpf-^A J Q^ (11) 

^3 

which now contains a Lagrange multiplier Xj for each of the locally conserved quan- 
tities. The existence of additional local conservation laws therefore severely restricts 
a possible thermalization of the system leading to additional constraints which are 
incorporated by the Lagrange multipliers in Eq. (11). Experimental indications for 
such constraints have been seen in realizations of the Lieb-Liniger model in ultracold 
gases. 16 



3. Low-energy description of the XXZ model 

In the following sections, we want to concentrate on one of the simplest integrable 
lattice models, the XXZ model 

H = J E "a ( C h+i + h - c ) h ( c h 2 ) + A ( Hj ~ 2 ) ( nj+1 2 



3 = 1 



(12) 

Here c (c^) annihilates (creates) a spinless fermion, J gives the energy scale, A char- 
acterizes the nearest neighbor density-density interaction with the density operator 
rij = CjCj. N is the number of sites and the boundary conditions might be either 
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periodic (sum runs up to N with c$ +1 = c^) or open (sum runs only up to TV — 1). 
h acts as a chemical potential. With the help of the Jordan- Wigner transformation 

S? -►„,.-!, S+^(-iy c]e^, SJ ^(-iy Cj e-^, (13) 
where 4>j = X^=i n l we can also express this model in terms of spin-1/2 operators 

N(N-l) 

H = J [S*S* +1 +S]Sy +1 +AS*S* +1 -hS*]. (14) 

3=1 

For both kinds of boundary conditions the XXZ model is integrable by Bcthe 
ansatz. 17,25-28 The Luttinger liquid approach is applicable in the critical regime 
which is given by —1 < A < 1 for h = 0. In general, the range of anisotropics 
for which the model is critical depends on the applied magnetic field h. In the free 
fermion case, the model (12) is easily solved by Fourier transform leading to 

H a = e p c! p Cp with e p — — J(cosj3 + h) . (15) 
p 

where we have set the lattice constant a = 1. The allowed momenta arc given by 
p = 2-Kn/N with n = 0, • • • , N — 1 for periodic boundary conditions (PBCs) or 
p = nn/(N + 1), n = 1, • • • , N for open boundary conditions (OBCs). In Sec. 4 we 
will briefly discuss the Bethe ansatz solution of this model for OBCs. 

Let us first revisit the derivation of an effective low-energy description, the Lut- 
tinger theory, by bosonization following Refs. 4,29,30. First, we replace the fermionic 
operators in the continuum limit by two fields ipR,L defined near the two Fermi 
points ±kp — ± arccos(— h): 

Cj -> V (x) = e ikFX i> R (x) + e- ikFX ip L (x) . (16) 

In a second step, we use standard Abelian bosonization to write the fermion fields 

as 

iP R . L (x)~ -^=e-^^ x \ (17) 

where a <~ kp 1 is a short distance cutoff. Instead of working with the left and right 
components 4>r,l we can define a bosonic field <j) and its dual field 9 by 

(18) 



y/2 ' V2 

which satisfy the standard bosonic commutation rule [<fi (x) , d x >0 {x')\ = iS (x — x'). 

If we bosonize the kinetic energy term of Eq. (12) keeping only the lowest order 
we obtain 

i-L ... i-L 



H^ n = iv F J dx (: tf R d x il> R : - : ^ L d x ^ L :) = ^ J dx [(d x <j> R ) 2 + (8 X 



(19) 
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where :: denotes normal ordering and vf — J sin hp is the Fermi velocity. This 
approximation corresponds to a linearization of the dispersion e p at the Fermi points 
±kp. In this case the bosonic model is quadratic in d x (f>^L- Corrections to the 
kinetic energy appear due to band curvature. Including these curvature terms, we 
can write the expansion of the dispersion e p near the two Fermi points as 



R,L 



, k 2 7 fc 3 



(20) 



where k = p =p kp for the right or left movers, respectively M = (J cos kp)" 1 
is the effective mass and 7 = Js'mkp. Note that the inverse mass M" 1 vanishes 
in the particle-hole symmetric case h = 0. In this case, the curvature correction 
is cubic in momentum. Bosonization of the /c 2 -term leads to a correction cubic 
in d x (j)B„L, whereas the term cubic im momentum gives a quartic correction in 
terms of the bosonic fields. Cubic and quartic terms in the bosonic operators will 
also arise from the interaction term in Eq. (12). The scaling dimension of these 
terms is 3 respectively 4 so that they are formally irrelevant. The interaction will, 
however, also produce additional marginal terms, quadratic in the bosonic fields, 
which together with (19) lead to the exactly solvable Luttinger model 



H 



LL 



dx 



{d x (j> R f + (d x (f> L ) 2 



92 
TTVf 



d x (j>Ld x (i>R \ . (21) 



Here 52 = 9a = 2JA[1 — cos(2fcp)] = 4JAsin 2 kp are interaction parameters. The 
Hamiltonian (21) can be rewritten in the form 



Hll = IJ dx vK (d x ef + ^ (d x £f 



(22) 



where v (the renormalized velocity) and K (the Luttinger parameter) are given by 

2A 



v = vp\ 



K 



} + 2-kvf) Gttuf) ~ VF ' ' 



sin kp 




(23) 
(24) 



Expressions (23) and (24) are approximations valid in the limit |A| <C 1. In Sec. 5 
we will review how these parameters in the Luttinger liquid Hamiltonian can be 
fixed exactly for arbitrary interaction strengths — 1 < A < 1 using the Bethe ansatz 
solution. 

The Luttinger parameter in the Hamiltonian (22) can be absorbed by performing 
a canonical transformation that rescales the fields in the form (f> -4- \[K(\> and 
6 -> 9/y/K leading to 



H L l = - I dx 



{d x 0f + (d x 0f 



(25) 
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We can also define the right and left components of these rescaled bosonic fields by 
These are related to 4>r,l by a Bogoliubov transformation. 



(26) 



3.1. Irrelevant operators in the finite field case 

The leading irrelevant operators stem from the fc 2 -term in Eq. (20) and give rise to 
dimension three operators ~ {d x 4>R,L,) 3 - Similar terms will also arise by bosonizing 
the interaction term. Instead of deriving these terms from the microscopic Hamil- 
tonian, we can introduce them phcnomenologically by considering the symmetries 
of the problem. In particular, the low-energy effective theory has to be symmetric 
under the parity transformation <pL — > <p R , <f> R — > <f>L, and x — > —x. We therefore 
can parametrize these terms as 30 

/2tt 



SH 



dx 



|?7_ (d x ip L f - {d x ip R y 



+V+ (dx^L) 2 d x ip R - {d x tp R f d x tp L ^ | 

(27) 

We will see in Sec. 5 that we can relate the amplitudes rj± to quantities which are 
known from the exact solution. The derivation of these terms starting from the 
microscopic Hamiltonian, on the other hand, would only allow us to obtain the 
coupling constants to first order in A with 30 



V- 



1 

M 



1 



2A 



sin ki 



V+ 



3A 
'ttM 



sinfc 



F- 



(28) 



From this expansion we see that (a) both terms vanish in the limit h — > where 
M _1 — > 0, and (b) that the term mixing right and left movers parametrized by rj + 
is only present in the interacting case, A^O. 



3.2. Irrelevant operators for zero field 

In the particle-hole symmetric case, h = 0, the first correction to the linear disper- 
sion relation is cubic in momentum, see Eq. (20). Instead of bosonizing this term 
starting from the microscopic Hamiltonian (12) we again introduce the correspond- 
ing terms in the bosonic model based on symmetry arguments. The dimension four 
operators allowed by symmetry can be parametrized as 



SH 



ttC- 
12 



: {d x <fR) 2 : : {d x f B f : + : (d x f L f : : (d x f L f 



+ 



■ (d x (pn) 



■ {d x (p L f ■ + 7i"C3 [: (d x (p R ) 3 : : d x <p L : + : (d x Lp L f : : d x f R :] 



(29) 



The explicit bosonization of the corresponding band curvature and interaction terms 
yields the coupling constants again only to lowest lowest order 



c_ » - J l + - 

7T 



AJ 



Cs = o. 



(30) 
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In addition, the Umklapp scattering term SHu ~ c 4lkFX '$> R (x) 1 & l(x)^ r (x + 
1)* l(x + 1) + h.c. is commensurate in this case, Akp = 2n, and therefore has 
to be kept in the low-energy effective theory. Bosonizing this term leads to 

5Uu = Xcos(AV^K4>) , (31) 

and to lowest order in A we have A = JA/(2tt 2 ). For OBC, there is also an irrelevant 
boundary operator allowed 

5H B ~ (5(x) + 5(L))(d x 0) 2 . (32) 

Finally, we want to consider the Luttinger model (25) with an additional small 
magnetic field Sh added, ignoring the irrelevant terms 

r-L 



By performing a shift in the boson field 



{d x <t>f + {d x 9f--\-8hd x(l 



(33) 



K X Sh (34) 

7T V 

we return to the quadratic Hamiltonian (25) with an additional constant shift 
—LK{8h) 2 /(2irv). The bulk susceptibility per site is therefore given by 

Xbuik = K/irv . (35) 

This result does not only hold for h = but also for any finite field ho at which we 
want to calculate \ with K — K(A, ho) and v — v(A, h ). 

4. The Bethe ansatz solution 

To exactly solve the interacting system for PBC or OBC, one can use the coordinate 
Bethe ansatz. 17 ' 25-28 Here we want to review very briefly some of the essential 
results needed to fix the parameters in the Luttinger model and refer the reader to 
Rcfs. 28,31-33 for a more detailed discussion. The coordinate Bethe ansatz starts 
from the fully polarized state ('the vacuum') and one derives coupled eigenvalue 
equations H\M) = E\M) for states \M) with M spins flipped. The eigenenergies 
can then be written as 

M ( N - 1 \ 
E = jJ2cosk j + JA{— M\ . (36) 

j=i \ / 

The structure is similar to the non-interacting case, however, the momenta kj are 
shifted from their positions for A = 0. They can be determined from a set of coupled 
nonlinear equations. In the thermodynamic limit, a single integral equation for the 
density of roots p(x) is obtained which parametrizes the allowed momenta 

1 f B 

#(x, 7) + [#( x > 7) + 7T - 7) + !?(*, 2 7 )] = p{x) + J #{x- y, 2j)p(y) dy, 

(37) 
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where 

#(x^) = -—^ (38 ) 

7r cosh zx — cos 7 

and we have set A = cos 7. Eq. (37) is the integral equation for the XX Z chain with 
OBC in the thermodynamic limit, including the boundary correction ~ 0(1/N). 
Omitting the 1/N correction this is the standard integral equation for PBC. 28 The 
integral equation contains an unknown parameter B and an unknown function p(x). 
It can be solved analytically by Fourier transform for B — 00 and one finds (ignoring 
the 1/N correction) 

l y ' 27Cosh7ra;/7 v ; 

The magnetization per site m and the ground state energy per site e for general 
B are given by 



to = 1/2-/ p(x)dx + 1/(2N) 

J-B 

, , Jsin7 f B , , J f 2-cos7\ 
e = -hs* ^ J ${x, 1 )p{x)dx + - (cos 7 + ^1 . (40) 

Inserting (39) into Eq. (40) one finds that B = 00 corresponds to to = 0, i.e., to 
the case of zero magnetic field h — O.ln general, B = B(h) and the dependence on 
magnetic field has to be determined numerically. This can be achieved by using the 
stationarity condition 

±[e(h)-e(h = 0)}=0. (41) 



5. Fixing parameters of the Luttinger model using integrability 

One of the main motivations to apply Luttinger liquid theory to integrable mod- 
els is that parameters in the Luttinger liquid theory such as velocities, Luttinger 
parameters, coupling constants of irrelevant operators and prefactors of correla- 
tion functions which usually are non-universal and therefore unknown, can often 
be determined in the case of an integrable model. 31, 34-37 This makes it possible to 
obtain parameter-free results at low energies. The general idea is to calculate static 
observables at zero or finite temperatures exactly using the Bethe ansatz and to 
compare with results obtained within Luttinger theory. In the following we briefly 
review this method to obtain the velocity and Luttinger liquid parameter, Sec. 5.1, 
and the coupling constants for band curvature and Umklapp terms for the XXZ 
model, Sec. 5.2. 



5.1. Velocity and Luttinger liquid parameter 

To obtain the velocity of elementary excitations, we have to consider the change in 
energy when replacing the ground state distribution of roots by a distribution which 
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contains an excitation near the Fermi points. In the case of finite magnetic field, the 
obtained Bcthc ansatz equations can only be solved numerically Here we want to 
restrict ourselves to the zero field case. Replacing the ground state distribution (39) 
in the expression for the energy (40) by the distribution containing an excitation 
gives the energy in terms of the momentum change. 28 This allows to read off the 
spin velocity 

dE 



dk 



J7rsin7 Jn \f \ — A 2 
k=k 2 7 2 arccosA ' 



The spin velocity therefore increases from vp = J (remember that we have set a = 1) 
at the free fermion point to v = Jtt/2 at the isotropic antiferromagnctic Hcisenberg 
point. Conversely, the velocity vanishes, as expected, for A — > — 1 corresponding to 
the isotropic fcrromagnet. 

To determine the Luttinger parameter K it is easiest to calculate the bulk sus- 
ceptibility Xbuik(^) = dm/dh using Eq. (40). To do so, B(h) is required. For finite 
magnetic fields this again requires a numerical solution. For infinitesimal fields, on 
the other hand, B(h) can be determined analytically 28,32 ' 33 and 

1 7 1 

Xbulk = ~f 7 \ — : = TT! T ( 43 ) 

J (tt — j)tt sin 7 2v(n — 7) 

is obtained. Comparing with Eq. (35) we find 

K = - = - . (44) 

2(77-7) 2(7T - arccos A) 

Therefore K = 1 for the free fermion model and K = 1/2 at the isotropic antifer- 
romagnetic Heisenberg point. Eq. (42) and Eq. (44) agree to first order with the 
expressions (24). The velocity and Luttinger parameter, both for zero and finite 
fields, as obtained from the Bethe ansatz solution, are shown in Fig. 1. 



5.2. Coupling constants of irrelevant operators 

Next, we want to review how the coupling constants C± i3 , Eq. (29), and A, Eq. (31), 
can be fixed in the zero field case, and the coupling constants rj±, Eq. (27), for finite 
field. The zero field case has been first considered by Lukyanov 31 and analytical 
formulas for the coupling constants have been obtained. The finite field case has 
been treated in Refs. 30, 38 leading to formulas which require a numerical solution 
of the Bcthc ansatz equations. 



5.2.1. The zero field case 

The simplest way to determine the Umklapp scattering amplitude A is to consider 
an open chain with a small magnetic field added. 33 In the low-energy description 
this means that we have to consider (33) with the Umklapp term (31) added. We 
can then again perform the shift (34) . This brings us back to the standard Luttinger 
liquid Hamiltonian (25) and the magnetic field now appears in the Umklapp term 
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E (l) _ 



Fig. 1. (a) Velocity v/ J and Luttinger parameter K for magnetization m = as a function of A. 
Velocity (b) and Luttinger parameter (c) as a function of magnetization m for A = 0.2, 0.5, 0.7, 0.9 
(in arrow direction). 

(31). In first order perturbation theory in Umklapp scattering we then find the 
following boundary correction to the ground state energy 39 

\J^dx(cos(^V^K(b+ 41 ^-^ . (45) 

Here (• • • )o denotes the correlation function calculated for the free boson model. 
For PBC this correlation function would vanish, however, for OBC we obtain 

cos 

Note that this is a 1/N correction to the ground state energy per site e = E/N. By 
partial integration we can split of the convergent part and find 

4K-1 



E (i,canv) = _x(2K)* K r(-4K)sm(2Kn) I - I . (47) 



At the same time, we can apply the Bethe ansatz to analytically calculate the 
so-called boundary susceptibility xb given by \ = xtmik + Xb/N + 0(N 2 ) to lead- 
ing orders in /i. 32 ' 33 The amplitude A of the Umklapp term can now be found by 
comparing the exact result for xb with Eq. (47). This leads to 

4K-2 



_ 2KT(2K)sinir/2K 
7rr(2 - 2K) 



n 1 

2v^r(i 



4K-2 
i K 



2K-1 



(48) 



In Ref. 31 this result has been obtained first by calculating the bulk correction to 
the ground state energy. Note, however, that this requires second order perturbation 
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theory in the Umklapp scattering. Particular care has to be taken when considering 
the isotropic antiferromagnet, A = 1. In this case, Umklapp scattering becomes 
marginally irrelevant and A has to be replaced by a running coupling constant 
which depends on the length scale the system is considered at. In general, both 
the length of the system and temperature will be of importance and the running 
coupling constant g(L,v/T) can be introduced by the replacements K — > (1 + g)/2 
and A — > — g/A. An explicit solution of the renormalization group equations for g 
is only possible if one of those two length scales dominates. In the thermodynamic 
limit, for example, this scale will be set by temperature alone and one finds 31 

l/ fl + ln(5)/2 = ln(T /T) (49) 

with T = y / 7r/2e 1 / 4+ ^ where 7 is the Euler constant. The scale T has again been 
fixed by comparing with the Bethe ansatz result for the bulk susceptibility Xbuik in 
the isotropic case. 

Here integrability has been used to fix a coupling constant. The conservation 
laws underlying integrability discussed in Sec. 2 can, however, have an even more 
profound effect. 13,14 For the XX Z model the first of the non-trivial conserved 
quantities is the energy current J E = Q3 given by 

jE 7-2 \ " f ay rrz ax ax nz ay ■ a / ax ay az az ay ax \ 

j 

+A(S*_ 1 S*S? +1 - S]_^S^ +1 )] . (50) 

The latter is defined by the continuity equation of the energy density at zero field 

jf + i-jf = -dtH i =i[H j ,lT\, (51) 

where H = £ i Hj is the Hamiltonian (14) with h = 0, PBC, and J E = Ejif- 
The energy current operator for the Luttinger model can be obtained from (51) by 
taking the continuum limit. This leads to 

Jo = J dxj^{x) = V — J dx {d x (p R ) 2 - {d x (p L f = -v 2 J dxd x (f>d x 9. (52) 

This operator is conserved, i.e., [J E , Hll] = 0. The irrelevant operators (29) lead to 
a correction of the energy current which can again be calculated using the continuity 
equation (51). To first order one finds 



S J = ttv I dx 



3 



(d x (pR) 4 - (d x ip L ) 4 +2( 3 (d x (p R f d x tp L - (d x tp L f d x (p R j 



(53) 

For J E = J E + 5J E to be conserved as required by integrability, we have to require 
that [J E , H] = [J E +5J E , H LL +5H] = up to the considered order. Since [J%, H] = 
[Jq,H ll + SH] = this implies that [5J E ,H LL ] = 0. The (--term in Eq. (53) 
docs not mix right and left movers and therefore obviously commutes with H LL . 
The (3-term, on the other hand, does mix the two modes and therefore does not 
commute with H LL . Integrability therefore implies that £ 3 = 0. We see that apart 
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from determining the precise values of coupling constants in the low-energy effective 
theory, there is a more fundamental consequence: Integrability corresponds to a fine 
tuning of the coupling constants such that the local conservation laws are fulfilled. 
In particular, terms which are in general allowed by symmetry might be absent. 

We are left with only two amplitudes, £±, for the dimension four operators. Let 
us briefly review how they can be fixed as well. Using Eq. (26) we can express both 
terms by the boson field (f> and the dual field 9. Now performing again the shift (34) 
for a small applied magnetic field Sh we find a first order correction to the ground 
state energy per site 

The ((5/i) 4 -term in the ground state energy can also be calculated analytically by 
Bethe ansatz. 31-33 One finds that the result consists of two distinct, additive, con- 
tributions. One of those vanishes at the free fermion point and is therefore associ- 
ated with the C + -term in the low-energy effective theory which mixes right and left 
movers. The other term then determines £„ leading to 31 

v r (4K-2) (4^-2) v ( ttK 



r ( 3 r 3 ( — ) 



4tvK r f _3 \ r3 ( k_ V 2tt \2K 



C+ = -^-tan — - . (55) 



The dependence on anisotropy of all three coupling constants A, £± is shown in 
Fig. 2. Note that the amplitude ( + diverges for K = (2n + l)/4n with n E N. 
Corrections to observables calculated in perturbation theory in the irrelevant oper- 
ators are, however, usually finite. What happens is that at these special points the 
scaling dimensions of different irrelevant operators coincide and the two diverging 
amplitudes 'conspire' to produce a finite result. This point has been discussed in 
detail in Ref. 33 using the susceptibility as an example. 

5.2.2. The finite field case 

For a finite magnetic field (B < 00), the Bethe ansatz integral equation cannot be 
solved analytically. However the amplitudes n± can be related to changes in the 
Luttinger parameter and velocity when changing the field. 30 ' 38 This allows for an 
accurate numerical determination of these parameters. 

The basic idea is again quite simple. We consider the Luttinger liquid Hamilto- 
nian at some finite magnetic field h . This means that our left and right modes live 
near Fermi points ±kp ^ ±7r/2. Now we apply an additional small magnetic field 
5h which we can take care of by the boson shift (34) . If we now calculate the free 
energy we obtain 

ttT 2 K{h) . t:u . 2 
f = -^)-2^ JSh) (56) 

=X(h)/2 
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Fig. 2. (a) as function of anisotropy A. The amplitude diverges whenever the argument of 
the tan-function in Eq. (55) is ±7r/2(mod 2w). (b) f_, and (c) A, Eq. (48) as a function of A. 



with h = ho + Sh. The interaction parameters K{h) and v(h) can be determined 
numerically as described in Sec. 4. Now we can expand (56) in Sh and obtain to 
lowest order 



Sf 



irT 2 dv 



6v 2 (h )dh 



Sh; 8x 



h—ho 



K(h ) 
irv(ho) 



1 dK 

Kim 



h—ho 



1 dv 
v dh 



h—ho 



Sh (57) 



These corrections have to stem from the dimension three operators (27). The second 
approach therefore is to keep K — K{h — h^), v = v(h = ho) fixed and to perform 
the shift (34) also in (27). Calculating again the free energy by standard techniques 
we now find 



Sf = (3r?_ 



7r \/KJShT 2 . 



-v+) 



K 3 / 2 (JShf 



K 3 / 2 J 

Sx = {7]+ - r] - ) mrW) 6h 
(58) 

A comparison of (57) and (58) yields two equations from which one obtains 30 ' 38 



6irv 3 (ho) 



dv 



dK 



3v 2 dK 



V ~ = K T i 2 dh + 2K 3 ! 2 ~dh' " 1+= 2KV 2 ~dh' [ > 

From this result a few general conclusions can be drawn. For a free model, rj + should 
vanish because a mixing of right and left movers is then impossible. This is indeed 
the case since K does not change when applying a field in this case. Furthermore, 
rj + will also be absent for models such as the Calogero-Sutherland model where 
K remains independent of the applied field even in the interacting case. The r}_- 
tcrm, on the other hand, is already present in a non-interacting system due to band 
curvature with r\_ = M^ 1 for A = 0. The result (59) of course also agrees with the 
expansion for small A, Eq. (28). 
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The parameters T)± are shown in Fig. 3 for different anisotropics A as a function 
of magnetization m. While rj± — > for m — > if K > 5/8 as expected from the weak 




-o. 
- + -o. 

-0. 
-0. 



1 1 1 1 1 1 ' 1 1 


: (b) 


i , i 


i 







2 
4 
6 
8 

-h.5 



-0.4 



Fig. 3. (a) r\— as a function of magnetization m for A = 0.2, 0.5, 0.7, 0.9 (in arrow direction), (b) 
Same for r;_|_. Note that r}± — > — oo for m — > in the case A = 0.9 leading to a sign change of rj—. 

coupling expansion (28) both parameters diverge, on the other hand, for K < 5/8 
(A > cos(7r/5) « 0.81). This behavior is discussed in more detail in Ref. 14 and we 
remind the reader that divergencies also occur in the amplitude £ + for h = 0. 

6. Applications 

We now want to consider a few examples where the low-energy effective theory has 
been used to obtain parameter-free results for several important observables. 

6.1. Impurities, Friedel oscillations and nuclear magnetic 
resonance 

One of the best known realizations of the spin- 1/2 antifcrromagnetic Heisenberg 
chain is the cuprate Sr 2 Cu03. 6 In this system excess oxygen dopes holes into the 
chain which seem to be basically immobile. 7,40 Effectively, this leads to randomly 
distributed non-magnetic impurities which cut the spin chain into finite segments. 
The magnetic properties are therefore determined by an ensemble of finite chains 
of random length N with OBC. 

In a chain with OBC, translational invariance is broken leading to a position 
dependent local susceptibility 

Xi - J-<^>fc=o = i<£*S£t>/,=o (60) 
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where T is the temperature and S^ ot — . In order to calculate Xj m t ne 

low-energy limit, we can express the spin operator in terms of the bosonic field 

-d x $ + c(-iy cos V4irK<i> . (61) 

Here the pref actor of the uniform part is fixed by the condition J^. Sj = S^ ot . 4,29 
The amplitude of the alternating part, on the other hand, can be fixed with the 
help of the Bethe ansatz solution. The techniques required are, however, much more 
involved than the ones reviewed in the previous section. In particular, one finds that 
c = ^A z /2 with A z as given in Eq. (4.3) of Ref. 34. 

Using Eq. (61) we can write Xj = X um + ( — tyxf- The uniform part x unl for 
the Luttinger model is given by 



d 2 f 

uni w J 

x ' dh 2 



7TV 



1 J2s z *$1 ex P [ 2KL ! 



h =o LT Es ,exp[-^S|] 

For LT/v — > and L even Xs ~ jzp cxp [— 2 klt \ wnereas f° r L odd Xs ~ (4LT) -1 . 
For LT/v — > oo the thermodynamic limit result (35) is recovered. Note that this 
zeroth order result is position independent and shows scaling with LT. Corrections 
to scaling occur due to the irrelevant bulk and boundary operators. For < A < 1 
the leading bulk irrelevant operator is due to Umklapp scattering (31). This leads 
to a first order correction in the free energy 

S h = lj o da:<cos(4V^A : ^)>s,exp(-87r«'<#)oBc) (63) 
where we have used the mode expansion for OBC 

/ 7T fit ' x sminnx/L) / ;„„vt ;„„*>t + \ . 

*0M) = + ^-S LtJ - + £ (e— -a„ + e—oi) (64) 

n— 1 v 

to split the expectation value of the Umklapp operator into an S z (zero mode) 
and an oscillator part. Furthermore, we have used the cumulant theorem for the 
oscillator part. It is now straightforward, although a bit tedious, to evaluate the 
two parts of (63). From this the correction to the uniform part of the susceptibility 
in first order in Umklapp scattering can readily be obtained. 11 
The boundary operator (32) yields a further correction 11 



s x r- nvb 



2KT 2 L 3 



(7TI.S 2 
V Sle-5KT* 



Q~ 2KLT 



(e ; 



(65) 



In the thermodynamic limit, Eq. (65) reduces to ^x™ 1 — > Kb/(-KvL). The field 
theory result in this limit can be compared with the calculation of the boundary 
susceptibility based on the Bethe ansatz 32 ' 33 and the proportionality constant b can 
be fixed 

fe = 2" 1/2 sin [wK/(4K -2)] /cos [tt/(8K-4)} . (66) 
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To first order in Umklapp scattering and in the dimension three boundary operator 
a parameter-free result for \ nnl can therefore be obtained. 

The alternating part of the susceptibility (60) can be written as xf = 
7^ exp(— 2irK(<fi(j>) osc )(cos V AirK(f)) s z where we have again split the correlation func- 
tion into an oscillator and a zero mode part using the mode expansion. The calcu- 
lation is now completely analogous to the calculation of the correction (63) leading 
to 



Xf = 



N + l 



K 



rf K (e 



£ m msin[27rmj/(7V + l)]e- 



-Kvm 2 /(2KLT) 



7 1 



-■Kvm 2 /(2KLT) 



(67) 

Here rj(q) is the Dedekind eta- function, and 9i(u, q) the elliptic theta function of the 
first kind. In the thermodynamic limit, L — Na oo, where we have reintroduced 
the lattice constant a for clarity, we can simplify our result and obtain 

2cK x 



X, 



A" 



with x = ja. This agrees for the isotropic Heisenberg case, K = 1/2, with the result 



in Ref. 41. In Fig. 4 the parameter-free formula for \j 
to Quantum-Monte-Carlo data. 12 



X" 



' + (— l) J xf i s compared 




Fig. 4. Local susceptibility for a finite open XX Z chain with length N, A = 0.3 at temperature 
T/J = 0.02. The squares denote the result of the parameter-free field theory formula, the circles 
are results obtained by Quantum Monte Carlo calculations. 12 



The position dependent susceptibility is directly measured as Knight shift in 
NMR. The hyperfine interaction couples nuclear and electron spins and the Knight 
shift of the nuclear resonance frequency for a chain segment of length N is given by 
Kj N ^ — {le/ln) J2j> Xy \ where 7 e (7„) is the electron (nuclear) gyromag- 
netic ratio, respectively. The hyperfine interaction is short ranged so that usually 
only A and A ±x matter. For a random distribution of non-magnetic impurities 
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within a chain the NMR spectrum reflects the distribution of Knight shifts for 
an ensemble of spin chains with random lengths. This leads to rather complicated 
NMR spectra 7 ' 40 whose properties can be fully understood using the parameter- free 
results for the susceptibility discussed above. 12 



6.2. The spin-lattice relaxation rate and transport 

The spin current (or particle current in the fermionic language) for the XXZ model 
is defined by 

J = -yE (S+Si +1 S+ +1 Sf) « J dxdj (69) 

Whether or not the integrable XXZ model supports ballistic transport at finite 
temperatures has been the topic of a long-standing debate. 13 ' 14 ' 20 ' 42-44 As discussed 
in Sec. 2.1 ballistic transport is signalled by a non-zero Drude weight and related 
to the part of the current which cannot decay due to conservation laws, see Eq. (7). 
For finite magnetic field, the Mazur inequality indeed immediately yields a non-zero 
Drude weight. 20 In this case the conserved energy current (52) becomes 

Jo 



v 2 J dx d x 6 d x <j> - hv^j^- J dx d x 9 = J Q E + hj (70) 



where we have used again the shift in the boson field, Eq. (34). The equal time 
correlations in (7) can now be evaluated for the Luttinger model (25) and 14 

For T//i4 the Mazur bound obtained from the overlap with J E saturates the 
exact zero temperature Drude weight D{T = 0) = vK/2ir. 45 Furthermore, one 
can also use the Bethe ansatz to calculate the Mazur bound -DMazur exactly. The 
obtained result agrees with (71) up to temperatures of order J. 14 

For zero magnetic field, however, the overlap between all local conserved quan- 
tities Q n of the XXZ model which can be constructed from the transfer matrix (4) 
and the current J vanishes, because the Q n are even under particle-hole transfor- 
mations while J is odd. Recently, a quantity — not related to the conserved quan- 
tities obtained from the Bethe ansatz solution — has been constructed for an open 
XXZ chain which is conserved up to boundary terms. 46 This quantity seems to 
protect part of the current in the thermodynamic limit, a view which appears to be 
supported by new numerical data. 47 As in the finite field case, the correct picture 
therefore seems to be that at finite temperatures a diffusive and a ballistic trans- 
port channel coexist. 13 For temperatures T/J e [0.2,0.5] and h = 0, the Drude 
weight, however, seems to be much suppressed compared to its zero temperature 
value known exactly from Bethe ansatz. Furthermore, the Drude weight seems to 
vanish completely at finite temperatures in the isotropic case, A = 1. We therefore 
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ignore a protected part of the current for now and first concentrate on the diffu- 
sive channel. The corrections due to a possible conserved part of the current are 
discussed at the end of this section. 

From an experimental point of view, it is of great interest to understand the 
intrinsic mechanism which leads to spin diffusion in the XX Z model. Spin diffusion 
has directly been observed in the spin lattice relaxation rate 1/Ti measured in NMR 
experiments as a magnetic field dependence 1/Ti ~ 1/v^- 8 ' 48 Both the spin lattice 
relaxation rate and the conductivity can be calculated within the Luttinger model 
from the retarded boson propagator 

For n ret (g,u;) = this is just the free boson propagator. In the zero field case, the 
leading irrelevant operator is the Umklapp term (31) and we will concentrate here 
on calculating the self-energy IT et (g, u) in first order in this perturbation. Further 
contributions to the self-energy will stem from the band curvature terms and are 
discussed in Ref. 14. Note, however, that only Umklapp scattering can give the 
bosons a finite lifetime and thus lead to diffusive transport. The calculation of the 
correlation function (72) in first order in Umklapp scattering is straightforward. 49 
We find 13 ' 14 

n" Ct (q,Lj) RJ -2i"fLJ (73) 

where the decay rate 7 is given by 

2 7 - SttKX 2 sm(4nK) ( —) § V(l/2 - 2K)T(2K) 1 ~ ^ cot(2^)T 8JY - 3 

\ v J v 7r 2 4A+i 

(74) 

in the anisotropic case and by 

2 7 = irg 2 T (75) 

in the isotropic case, A = 1. The running coupling constant g — g(T) is determined 
by (49). The Kubo formula directly relates the conductivity a(q, u) to the calculated 
bosonic Green's function 

a( 9 ,w) = -iw<#) ret (g,a;) (76) 

and due to the finite relaxation rate 7 at finite temperatures one find a Lorentzian 
for the real part of the conductivity 

<77) 

At the same time, also the spin lattice relaxation rate can be expressed by the same 
bosonic Green's function 

1 ^ 2KT I fdq |j%)| 2 Im {#>ret(g)We) . (78) 

±1 TTUJ e J ZTT 
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Here to e = ^sh is the electron magnetic resonance frequency. 14 If the hypcrfme 
coupling form factor A(q) picks out the q <~ contributions of the integral (78) as in 
the oxygen NMR experiment on Sr2Cu03 in Ref. 8 then 1/T\ and the conductivity 
<j'(q ~ 0, u>) are directly related. The experimentally considered spin chains are 
almost isotropic, Aal. Using the parameter- free result (75) we obtain the diffusive 
behavior 




T/\vl(J/T) 



(79) 



The only free parameters remaining depend on microscopic details of the considered 
compound. Both the exchange constant J and the hypcrfme coupling constant A{q) 
can be fixed by analyzing the susceptibility and performing a K — \ analysis respec- 
tively. 6 ' 8,40 Integrability therefore makes it possible to obtain an analytical result 
for the spin-lattice relaxation rate which includes the intrinsic relaxation processes. 
The comparison of the result (79) with experiment as shown in Fig. 5 demonstrates 
that Umklapp scattering seems to be the dominant source for relaxation and that 
other contributions, e.g., due to electron-phonon scattering are small. 




100 



200 300 400 
Temperature (K) 



500 



600 



Fig. 5. Spin-lattice relaxation rate 1/Ti for Sr2CuC>3. Experimental data (symbols) are compared 
to the field theoretical result. The only free parameters are the exchange constant J ~ 1800 K 
and the hypcrfme coupling tensor which have been determined experimentally. A more detailed 
analysis is presented in Ref. 14. 



Finally, we have to discuss the possible conservation of part of the current even at 
zero magnetic field. The quantity constructed in Ref. 46 contains terms active over 
several lattice sites for < A < 1. It is therefore non-trivial to bosonize this term 
and to consider the consequences of its (almost) conservation for the low-energy 
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effective theory in a similar spirit to our discussion of J E in Sec. 5.2. An alternative 
method discussed in Refs. 13,14,50 is to use a memory matrix formalism to calculate 
the self-energy (72) in the presence of conservation laws. Instead of Eq. (73) the 
self-energy then reads 

n «' = ^Tr|^ " ith '- OT <£-Vq>- <80) 

where Q is the conserved quantity with (JQ) ^ 0. As a consequence, the current- 
current correlation function at large times is now given by 

1 KvT 

l {J{t)J) = 4IT^) [v + ex ^- 2 ^ + f )*)] ( 81 ) 

where the first term is proportional to the Drude weight and the second term de- 
scribes the diffusive part. For A = 1 the Drude weight seems to be zero at finite 
temperatures 46 so that our theoretical analysis of 1 jT\ is not affected. In any case, 
the spin chains in Sr 2 Cu03 do not represent an integrable system and there is there- 
fore certainly no ballistic channel. It is known, for example, that in this compound 
also a weak next-nearest neighbor coupling J 2 ~ 0.1J exists. While this coupling 
destroys integrability it will only lead to a weak renormalization of the Umklapp 
amplitude A. 29 It is thus legitimate to use the results for the integrable model to 
analyze the experimental data for 1 /Ti . 

7. Conclusions 

Integrable gapless one-dimensional quantum models allow to test many predictions 
of Luttinger liquid theory. As such they have been vital in confirming the universal 
applicability of the latter. On the other hand, Luttinger liquid theory has also helped 
to understand integrable quantum models better. Except for the simplest integrable 
systems, such as free bosonic or fermionic particles and to some extent also Calogero- 
Sutherland type models, an exact calculation of correlation functions for arbitrary 
distances and times has not been achieved yet. A combination of Luttinger liquid 
theory and integrability then often allows to obtain a more complete picture. For 
Bethe ansatz integrable systems such as the Lieb-Liniger Bose gas, the (anisotropic) 
Heisenberg, the Hubbard, or the supersymmetric t — J model it is possible to fix 
the velocity of the collective excitations v and the Luttinger parameter K as a 
function of density and interaction strength. For the Luttinger model, (dynamical) 
correlations can then be easily calculated. 

Often one can go even one step further and determine the amplitudes of leading 
irrelevant operators acting as corrections to the Luttinger Hamiltonian as well as 
amplitudes of correlation functions, exactly. This article is by no means a complete 
review of all the results which have been obtained in this field. Instead, I have used 
the XXZ model as an example and have summarized how v, K as well as the 
amplitudes of the leading irrelevant band curvature terms and Umklapp scattering 
can be determined. As an application, I have shown that this allows to derive 
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parameter-free results for the local susceptibility in open Heisenberg chains and 
quantitatively explains Knight shift spectra which have been investigated by nuclear 
magnetic resonance for compounds such as Sr2Cu03. As a second application, I 
have presented results for the relaxation rate of the Luttinger liquid bosons at 
finite temperatures in first order in Umklapp scattering. From this a parameter- 
free formula for the spin-lattice relaxation rate 1/Ti can be obtained which is in 
excellent agreement with experiment. The same bosonic correlation function and 
the same relaxation rate determine, on the other hand, also the particle transport. 
In general, the XX Z model has coexisting ballistic and diffusive transport channels. 
The relative weight of each of these channels can be determined by combining the 
Luttinger model — keeping the dangerously irrelevant Umklapp term with known 
amplitude — with a memory-matrix calculation. 

There are many more interesting results which have not been covered in this 
review. One of the perhaps most fascinating recent developments is the so-called 
non-linear Luttinger liquid theory which allows to calculate dynamic response func- 
tions while taking band curvature into account. 51 This has lead to the discovery 
of new power laws near edge singularities. For integrable models the exponents of 
these power laws can be determined exactly. 52 A comprehensive overview about 
these developments has been given in a recent excellent review. 53 
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